System and method for shear wave elastography by transmitting ultrasound with subgroups of ultrasound transducer elements

ABSTRACT

Systems and methods for performing shear wave elastography using push and/or detection ultrasound beams that are generated by subsets of the available number of transducer elements in an ultrasound transducer. These techniques provide several advantages over currently available approaches to shear wave elastography, including the ability to use a standard, low frame rate ultrasound imaging system and the ability to measure shear wave speed throughout the entire field-of-view rather than only those regions where the push beams are not generated.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application is a continuation of U.S. patent application Ser. No. 16/023,782 filed Jun. 29, 2018, which is a divisional of U.S. patent application Ser. No. 14/397,395 filed Oct. 27, 2014, now U.S. Pat. No. 10,624,609, issued Apr. 21, 2020, which is the US National Phase of PCT/US13/63631 filed Oct. 7, 2013, which claims priority to U.S. Patent Application 61/710,744 filed Oct. 7, 2012, the entirety of all of which are incorporated by reference herein for all purposes.

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH

This invention was made with government support under EB002167 and DK082408 awarded by the National Institutes of Health. The government has certain rights in the invention.

BACKGROUND OF THE INVENTION

The field of the invention is systems and methods for ultrasound imaging. More particularly, the invention relates to systems and methods for shear wave elastography using ultrasound.

Shear waves can be used to evaluate the viscoelastic properties of tissue, which are sensitive biomarkers of tissue pathology. To perform two-dimensional ultrasound shear wave elastography, an ultrasound imaging system with a frame rate of several kilohertz is required to capture the fast moving shear waves. This is not feasible with conventional ultrasound imaging systems, which typically have a frame rate of less than 100 Hz and where two-dimensional images are formed line by line. Thus, it would be desirable to provide a system and method capable of performing two-dimensional shear wave elastography with a conventional ultrasound imaging system.

SUMMARY OF THE INVENTION

The present invention overcomes the aforementioned drawbacks by providing a method for measuring a mechanical property of an object using an ultrasound system having an ultrasound transducer that includes a plurality of transducer elements. At least one shear wave is induced in the object, and elastography data is obtained from the object. The elastography data is obtained by dividing the transducer elements in the ultrasound transducer into a plurality of subgroups of transducer elements; successively transmitting a set of ultrasound beams using a different subgroup of transducer elements, thereby successively transmitting ultrasound to different regions in the object; and repeating that process at a pulse repetition frequency such that each set of ultrasound beams is effectively repeated at the pulse repetition frequency. From the obtained elastography data, a mechanical property of the object can then be calculated.

It is another aspect of the invention to provide a method for measuring a mechanical property of an object using an ultrasound system having an ultrasound transducer that includes a plurality of transducer elements. At least one shear wave is induced in the object by dividing the transducer elements in the ultrasound transducer into a plurality of subgroups of transducer elements and transmitting a plurality of focused ultrasound beams, each focused ultrasound beam being transmitted by a single subgroup of transducer elements. Elastography data is obtained from the object by transmitting a detection ultrasound beam into the object and receiving echo signals in response thereto. A mechanical property of the object can then be calculated using the obtained elastography data.

The foregoing and other aspects and advantages of the invention will appear from the following description. In the description, reference is made to the accompanying drawings which form a part hereof, and in which there is shown by way of illustration a preferred embodiment of the invention. Such embodiment does not necessarily represent the full scope of the invention, however, and reference is made therefore to the claims and herein for interpreting the scope of the invention.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is an illustration of a plurality of ultrasound beam sets, each containing a plurality of ultrasound beams, being transmitted by an ultrasound transducer;

FIG. 2 is an example of a pulse timing diagram using the plurality of ultrasound beam sets of FIG. 1 ;

FIG. 3 is an example of aligning the time grids of different ultrasound beam sets using interpolation;

FIG. 4 is an example of a delay time in a shear wave front;

FIG. 5 is an example of another delay time in a shear wave front;

FIG. 6 is an illustration of a plurality of focused ultrasound push beams being simultaneously transmitted in a comb pattern;

FIG. 7 is an illustration of a focused ultrasound push beam that is rapidly translated along a lateral direction to generate shear waves; and

FIG. 8 is an illustration of a plurality of focused ultrasound push beams being simultaneously transmitted, in which each of the ultrasound push beams is generated with different sized subgroups of transducer elements.

DETAILED DESCRIPTION OF THE INVENTION

Described here are systems and methods for performing shear wave elastography using push and/or detection ultrasound beams that are generated by subsets of the available number of transducer elements in an ultrasound transducer. These techniques provide several advantages over currently available approaches to shear wave elastography, including the ability to use a standard, low frame rate ultrasound imaging system and the ability to measure shear wave speed throughout the entire field-of-view rather than only those regions where the push beams are not generated.

High-end, conventional ultrasound imaging systems can image several lines from a single transmission using parallel beam forming. The pulse repetition frequency (“PRF”) and pulse repetition interval (“PRI”) for ultrasound imaging are determined by the maximum depth of the image. Assuming an imaging depth of 30 mm, the round-trip ultrasound travel time is

$\begin{matrix} {\frac{2 \cdot {30\lbrack{mm}\rbrack}}{1.54\left\lbrack \frac{mm}{µs} \right\rbrack} = {{39\lbrack{µs}\rbrack}.}} & (1) \end{matrix}$

Adding beam forming and other overhead time delay, a PRI of 50 μs (PRF=20 kHz) is feasible. The inventors have discovered that for two-dimensional shear wave elastography, the detection PRF at each spatial location can be less than 1 kHz. That is to say, using 1 kHz as an example, it is sufficient to perform shear wave elastography detection at the same spatial location only once every 1000 microseconds. Therefore, shear wave elastography can be performed using parallel beam forming detection at multiple locations sequentially while maintaining a PRF of 1 kHz at each location.

Referring now to FIG. 1 , an example of a sequence of ultrasound beams generated by an ultrasound system to produce a desired PRF at each spatial location is illustrated. For example, the sequence of ultrasound beams includes N sets, S₁, . . . , S_(N) 102 each of M beams, B₁, . . . , B_(M) 104 formed in parallel. This sequence is then repeated at the pulse repetition frequency for a number, P, of desired repetitions.

By way of example, each set, S_(n), of ultrasound beams may include a plurality of beams, B_(m), that are formed in parallel. For example, each set of ultrasound beams may include four beams formed in parallel and twenty beam sets may be used for a total sampling of eighty different locations during each pulse sequence. Using this approach, shear wave motion can be detected within a two-dimensional region covered by the M·N spatial locations with an effective PRF of 1 kHz at each spatial location. Assuming a line spacing of 0.3 mm, which is typical for a 5 MHz linear array ultrasound transducer, two-dimensional shear wave elastography can be performed within an area that is 30 mm deep and 24 mm wide using this technique. In other implementations, each set of ultrasound beams may include only a single ultrasound beam. As a result, the region-of-interest in which mechanical properties can be measured will be smaller than with parallel beam forming.

An example pulse sequence timing for P repetitions of the pulse sequence illustrated in FIG. 1 is illustrated in FIG. 2 . At time t=0 μs, transmission of beam set S_(1,1) is used to track shear wave motion at beam locations defined by the beams, B_(m), that form beam set S_(1,1). For example, when beam set S_(1,1) is composed of four beams, shear wave motion at four locations will be tracked. At a preset time after beam set S_(1,1) is generated, beam set S_(2,1) is generated to detect shear wave motion occurring at beam locations defined by the beams, B_(m), that form beam set S_(2,1). For example, the preset time may be 50 microseconds. In the next cycle of the transmission, the first beam set is again transmitted. This beam set, S_(1,2), is therefore transmitted at time t=1000 μs for a PRF of 1 kHz. The pulse sequence is repeated until the P^(th) repetition is performed at time t=P·1000 μs.

The size of the two-dimensional shear wave elastography region is determined by the parallel beam forming capability of the ultrasound system and the required PRF at each spatial location. For example, six-beam parallel detection can allow a 36 mm image width if everything else does not change. It is contemplated that the required PRF may be higher in stiffer tissues. In this instance, the two-dimensional shear wave elastography region will be reduced in size in order to sustain the higher PRF. Therefore, implementations on traditional imagers may end up with a region-of-interest whose size changes with different applications. This result is similar to current two-dimensional color imaging methods in ultrasound scanners.

It is noted that although the descriptions provided above are made with respect to a one-dimensional array transducer, the concepts are also readily applicable to transducers with higher dimension, such as a two-dimensional array transducers.

It is noted that even though the PRF at each location is 1 kHz, different locations are sampled at different time grids. Referring to FIG. 2 , the four beam locations in beam sets S_(1,p) are sampled at times t=0,1000, 2000, . . . , P·1000 μs

The beams in beam sets S_(2,p) are sampled at times, t=50,1050,2050, . . . , (P·1000)+50 μs.

This small time delay among locations needs to be accounted for in shear wave elastography reconstruction. Two example methods for accounting for this delay time are provided below.

One example method for compensating for the delay time noted above is to interpolate the time signal at each location from PRF to N·PRF (e.g., from 1 kHz to 20 kHz) and to use the same time grid for all locations. As shown in FIG. 3 , the squares represent the time instances of ultrasound detection at locations covered by beam sets S₁ to S_(N). Each beam is sampled at a PRF of 1 kHz, but there is a time shift of 50 μs between adjacent beams. With interpolation back to 20 kHz for each beam, the time grids for each beam are aligned and synchronized. The interpolated points are represented by filled circles in FIG. 3 .

Referring to FIG. 4 , the triangular waves in (a) and (b) represent the shear wave signals detected at two locations x and y after interpolation such that both signals start at the same time. Direction x→y is the direction of shear wave propagation. The time delay, Δt, of the shear wave between (a) and (b) can be calculated using time-to-peak, cross-correlation, or other methods. The shear wave speed can therefore be calculated by

$\begin{matrix} {{c_{s} = \frac{\left| {y - x} \right|}{\Delta t}};} & (2) \end{matrix}$

where |y−x| is the distance between location y and location x.

An example of another method for compensating for the delay time noted above is to account for the delay time when calculating the shear wave speed. Referring to FIG. 5 , the triangular waves in (a) and (b) represent the shear wave signals detected at two locations, x and y, where the direction x→y is the direction of shear wave propagation. The time delay, Δt, of the shear wave between (a) and (b) can be calculated using time-to-peak, cross-correlation, or other methods. As explained above, the shear wave at different locations are detected at different time grids. Assuming that shear wave detection at location x starts at time t=0 as shown in (a), and shear wave detection at location y starts at time t=δ as shown in (b), the actual delay, ΔT, between the shear wave at locations x and y is ΔT=(Δt+δ)  (3).

By way of example, if x is located at one of the beams in beam set S₁, and y is located at one of the beams in beam set S₃, then δ=(3−1)·Δt=2·50=100 μs  (4).

The shear wave speed, c_(s), can then be calculated by

$\begin{matrix} {{c_{s} = \frac{\left| {y - x} \right|}{\Delta T}};} & (5) \end{matrix}$

where |y−x| is again the distance between location y and location x.

It is noted that tracking beam sets S₁, . . . , S_(N) do not need to be in a spatially sequential order; instead, the tracking beam sets can be placed in different spatial orders. For example, tracking beams set S₁ can be placed on left end of the image, whereas tracking beam set S₂ can be placed at right end of the image, with odd-numbered beam sets following sequentially inward to the center of the image from beam set S₁ and even-numbered beam sets following sequentially inward to the center of the image from beam set S₂.

In addition, two or more beam sets can be transmitted simultaneously to reduce the number of transmit-detection events required to cover a two-dimensional region of desired size. This approach has the benefit of increasing the effective PRF of detection at each beam set location. For example, beam sets S₁ and S₁₀ can be transmitted and detected simultaneously, followed by beam sets S₂ and S₁₁ simultaneously, until all N beam sets are covered. This process can repeat again in time.

It is also noted that each of the beam sets will be transmitted by a group of transducer elements. Different beam sets may have same or different transducer elements. For example, in a linear array or curved array ultrasound transducer, elements 1 through 32 may be used for beam set S₁, and elements 8 through element 40 may be used for beam set S₂, and so on. In another example of a phase array transducer with 64 elements, the same 64 elements may be used to transmit all beam sets from S₁ to S_(N) by steering different beam sets to cover different regions.

The above teaching is for calculating the shear wave speed by measuring the time delay between shear waves detected at two or more locations along the shear wave propagation direction. The same concept can be applied for other methods of shear wave elasticity imaging, such as direct inversion and frequency dependent dispersion analysis. In addition, this method can be applied to shear waves produced by physiological motion, such as cardiac motion; mechanical vibration; and ultrasound radiation forces from single focused beams, single unfocused beams, multiple focused beams, or multiple unfocused beams.

Having described applications for transmitting ultrasound detection beams using subsets of the available transducer elements in an ultrasound transducer, techniques for using subsets of transducer elements to deliver ultrasound push beams are now discussed.

When an ultrasound push beam is transmitted into an object, shear waves are generated and propagate outward from the push beam in opposite directions. Consequently, shear waves are not generated in the push beam region, which means that shear wave speeds cannot be measured in the region where the push beam is generated. In addition, shear waves produced by an ultrasound push beam attenuate quickly over a short propagation distance. As a result, shear waves produced by a single push beam can only image a small region-of-interest; thus, in these instances, multiple push-detect acquisitions are required to piece together an image with large field-of-view (“FOV”).

As described in co-pending PCT Application Publication WO2012/116364, entitled “Ultrasound Vibrometry with Unfocused Ultrasound,” which is herein incorporated by reference in its entirety, a comb-shaped set of unfocused ultrasound beams can be used to provide a full FOV, two-dimensional shear wave speed map together with one rapid data acquisition. This method is referred to as comb-push ultrasound shear elastography (“CUSE”).

In CUSE, multiple unfocused ultrasound push beams are used to produce shear waves within a tissue for shear wave elasticity imaging. Only one subset of transducer elements is used for each push beam; thus, multiple subsets of elements can be used for different spatial locations to simultaneously transmit multiple push beams. In CUSE, shear waves produced by each push beam can be treated as an independent realization of a single push beam.

Using CUSE, shear waves from different push beams interfere with each other and eventually fill the entire field-of-view (“FOV”). To achieve robust shear wave speed estimation, a directional filter can be used to extract left-to-right (“LR”) propagating shear waves and right-to-left (“RL”) propagating shear waves from the interfering shear wave patterns. A time-of-flight based shear wave speed estimate method may be used to recover local shear wave speed at each pixel from both LR waves and RL waves. A final shear wave speed map may then be combined from the LR speed map and RL speed map. Because comb-push pulses produce shear wave motions with high amplitude at all image pixels, including at the push beam areas, both shear wave speed at the “source free” areas and shear wave speeds at the push beam areas can be recovered.

Thus, CUSE enables a full FOV two-dimensional reconstruction of a shear elasticity map with only one data acquisition. To improve acoustic radiation force penetration and generate stronger shear waves into deeper tissue (i.e., liver and kidney), the previous CUSE method can be modified using focused ultrasound push beams. Using focused ultrasound push beams facilitates the generation of strong shear waves at locations deep within tissues. This ability to generate strong shear waves at deep tissue locations can lead to higher SNR for shear wave elasticity imaging compared to CUSE with unfocused beams.

In one configuration, illustrated in FIG. 6 , the transducer elements are divided into a number of subgroups, such as four subgroups, that each simultaneously transmits a focused ultrasound beam. This technique is referred to as focused CUSE, or “F-CUSE.” As noted, in F-CUSE, the transducer elements are divided into a number, N, of subgroups, with each subgroup containing one or more transducer elements.

As one example of the F-CUSE technique, a 128 element ultrasound transducer can be divided into four subgroups of 32 transducer elements each. In F-CUSE, all subgroups transmit focused ultrasound beams simultaneously to form a comb-push pattern ultrasound field. In one example, the duration of the push pulse beams can be on the order of 600 μs.

In another configuration, the transducer elements are divided into a number of subgroups with overlapping elements, and the subgroup that is used to transmit an ultrasound beam is rapidly changed along the lateral direction. The result of this “marching” of the subgroup of transducer elements being energized is to provide a focused ultrasound push beam at successively different horizontal locations. This technique may be referred to as marching CUSE, or “M-CUSE.”

As one example of the M-CUSE technique, shown in FIG. 7 , a 128 element ultrasound transducer can be divided into four, overlapping subgroups of 64 elements each. The first subgroup then transmits a single, focused push beam at a first time, t1. The duration of this push beam can be shorter than would be used in F-CUSE or traditional CUSE if it is desirable to control tissue or transducer heating due to repeated transmission using the overlapping transducer elements. As an example, the push beam duration in M-CUSE can be on the order of 200 μs. The push beam duration is selected based on considerations for how much overlap there is between subgroups of the transducer elements. For instance, the duration can be selected such that any given transducer element in overlapping subgroups is not energized for a consecutive duration that may lead to overheating in that element.

After the first push beam is transmitted, a second push beam is transmitted at a second time, t2, using the second subgroup of elements. This continues for the third and fourth subgroups. There is typically a small duration between consecutive push beams. For example, there can be a 15 μs delay between consecutive push beams. In this short duration of time, shear waves will propagate in soft tissue only about 0.45 mm, which is about 1.5 times the size of an individual transducer element in the transducer. Thus, in general, the amount of shear wave propagation between successive push beams is negligible for all subgroups after transmitting all of the focused push beams.

It is noted that each push beam of the F-CUSE or M-CUSE techniques may have a different number of transmit elements and may be focused at different depths. Additionally, each of the push events for the M-CUSE technique may include using more than one push beams. In the example shown in FIG. 8 , push beams 1 (with less elements) and 2 (with more elements) are simultaneously transmitted at time t₁, followed by push beams 3 and 4 transmitted simultaneously at time t₂. This combination push will produce strong shear waves at all depths and all lateral positions for shear wave detection and processing. It is also noted that push beams in M-CUSE can have different push durations and can be transmitted in an arbitrary order (for example, in a different order 1→4→3→2 compared to the order 1→2→3→4 illustrated in FIG. 7 )

Similar to the original CUSE method, both F-CUSE and M-CUSE can generate comb-patterned ultrasound push beams that induce a complicated shear wave field with interferences. Directional filtering described by Manduca et al. in “Spatio-Temporal Directional Filtering for Improved Inversion of MR Elastography Images,” Medical Image Analysis, 2003; 7(4): 465-473, can thus be used to separate the shear waves into multiple directions without interference so that robust shear wave estimates can be achieved at each imaging pixel within the FOV.

After comb-push transmission, a plane wave imaging mode can be used with all transducer elements delivering ultrasound to detect the propagating shear waves. Alternatively, the detection scheme described above can also be used.

The present invention has been described in terms of one or more preferred embodiments, and it should be appreciated that many equivalents, alternatives, variations, and modifications, aside from those expressly stated, are possible and within the scope of the invention. 

The invention claimed is:
 1. A method for measuring a mechanical property of an object using an ultrasound system having an ultrasound transducer that includes a plurality of transducer elements, the steps of the method comprising: a) controlling the ultrasound transducer to energize a plurality of subgroups of transducer elements so as to transmit a corresponding plurality of different focused ultrasound beams each focused at one of a plurality of different lateral locations, wherein the plurality of different focused ultrasound beams extend outward from the ultrasound transducer in a comb-shaped pattern; b) obtaining elastography data representative of the multiple different shear waves by transmitting at least one detection ultrasound beam into the object and receiving echo signals in response thereto; c) generating processed elastography data by processing the obtained elastography data to separate the multiple different shear waves into multiple different directional components to reduce effects of interference between the multiple different shear waves; and d) calculating a mechanical property of the object using the processed elastography data.
 2. The method of claim 1, wherein: the plurality of subgroups of transducer elements comprises a first plurality of subgroups of transducer elements and a second plurality of subgroups of transducer elements, the comb-shaped pattern comprises a first comb-shaped pattern associated with the first plurality of subgroups of transducer elements and a second comb-shaped pattern associated with the second plurality of subgroups of transducer elements; and step a) comprises: controlling the ultrasound transducer to energize a first plurality of subgroups of transducer elements so as to transmit a corresponding first plurality of different focused ultrasound beams each focused at one of the plurality of different lateral locations, wherein the first plurality of different focused ultrasound beams extend outward from the ultrasound transducer in a first comb-shaped pattern; and controlling the ultrasound transducer to energize a second plurality of subgroups of transducer elements so as to transmit a corresponding second plurality of different focused ultrasound beams each focused at one of the plurality of different lateral locations, wherein the second plurality of different focused ultrasound beams extend outward from the ultrasound transducer in a second comb-shaped pattern.
 3. The method as recited in claim 2, wherein at least some of the first plurality of different focused ultrasound beams and second plurality of different focused ultrasound beams are generated by spatially adjacent subgroups of transducer elements.
 4. The method as recited in claim 2, wherein at least some of the first plurality of different focused ultrasound beams and second plurality of different focused ultrasound beams are generated by spatially overlapping subgroups of transducer elements.
 5. The method as recited in claim 1, wherein step b) includes: i) dividing the transducer elements in the ultrasound transducer into a plurality of subgroups of transducer elements; ii) successively transmitting a set of ultrasound beams using a different subgroup of transducer elements, thereby successively transmitting ultrasound to different regions in the object; iii) repeating step ii) at a pulse repetition frequency such that each set of ultrasound beams is effectively repeated at the pulse repetition frequency.
 6. The method as recited in claim 1, wherein step c) includes applying a directional filter to the elastography data.
 7. The method as recited in claim 1, wherein step a) includes transmitting at least some of the plurality of different focused ultrasound beams to focus at different depth locations than other ones of the plurality of different focused ultrasound beams.
 8. The method as recited in claim 1, wherein calculating the mechanical property in step d) includes using at least one of a time-to-peak method, a cross-correlation method, a direct inversion method, and a frequency dependent analysis method.
 9. The method as recited in claim 1, wherein at least some of the plurality of different focused ultrasound beams are transmitted for a longer duration than others of the plurality of different focused ultrasound beams.
 10. The method as recited in claim 1, wherein some of the plurality of subgroups of transducer elements include more transducer elements than others of the plurality of subgroups of transducer elements.
 11. The method of claim 1, wherein the plurality of different focused ultrasound beams are focused at lateral locations in a first region in the object and the elastography data are obtained from a second region that is different from the first region.
 12. The method of claim 11, wherein the second region comprises lateral locations that are adjacent the first region. 